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Abstract 

We analyze the dynamics of rotary biomotors within a simple nano-electromechanical model, 
consisting of a stator part and a ring-shaped rotor having twelve proton-binding sites. This model 
is closely related to the membrane-embedded Fq motor of adenosine triphosphate (ATP) synthase, 
which converts the energy of the transmembrane electrochemical gradient of protons into mechan- 
ical motion of the rotor. It is shown that the Coulomb coupling between the negative charge of the 
empty rotor site and the positive stator charge, located near the periplasmic proton-conducting 
channel (proton source), plays a dominant role in the torquc-gcncrating process. When approach- 
ing the source outlet, the rotor site has a proton energy level higher than the energy level of the site, 
located near the cytoplasmic channel (proton drain). In the first stage of this torque-generating 
process, the energy of the electrochemical potential is converted into potential energy of the proton- 
binding sites on the rotor. Afterwards, the tangential component of the Coulomb force produces 
a mechanical torque. We demonstrate that, at low temperatures, the loaded motor works in the 
shuttling regime where the energy of the electrochemical potential is consumed without producing 
any unidirectional rotation. The motor switches to the torque-generating regime at high temper- 
atures, when the Brownian ratchet mechanism turns on. In the presence of a significant external 
torque, created by ATP hydrolysis, the system operates as a proton pump, which translocates pro- 
tons against the transmembrane potential gradient. Here we focus on the Fq motor, even though 
our analysis is applicable to the bacterial flagellar motor. 

PACS numbers: 87.16.A-, 87.16.Uv, 85.85.+j 



I. INTRODUCTION 
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Biological rotary motors, such as the Fq motor of ATP (adenosine triphosphate) synthase 
and the bacterial flagellar motor (BFM), convert energy of the transmembrane electrochem- 
ical potential into mechanical motion [l], I^. The gradient of the electrochemical potential in 
the living cells is maintained by the metabolic mechanism that translocates ions (here pro- 
tons) from the negative (or cytoplasmic) side of the membrane to its positive (periplasmic) 
side jl, y, y, 5|. A mechanical torque is generated when protons flow down the e 
ical gradient across the mitochondrial or cell membranes m, mi, 1121, 

n 

Bacteria employs this torque directly for chemotaxis [16], whereas in ATP synthase the ro- 
tational energy is transmitted to the Fi motor mechanically linked to the Fq portion. The 
out-of-membrane component Fi contains three sites, which harness the energy of mechanical 
rotations to catalyze the synthesis of ATP molecules. In the following, we concentrate on 
the Fq motor, although our theoretical analysis can be applied to the BFM as well, since 
these motors have many common features. 

A. Structure of the biomotor 



The Fq motor is embedded into a plasma membrane of bacteria or into an inner mi- 
tochondrial membrane. It consists of two components: (i) a stator or subunit a, tightly 
attached to the membrane, and (ii) a cylinder- shaped rotor (subunit cV which can freely 
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ISj the rotor 



rotate around its axis (Fig. 1). In the majority of models Baa 
contains several (10 to 15) proton-binding sites located in the middle of the membrane. The 
stator has two non-coUinear half-channels, which are perpendicular to the plane of Fig. 1. 
Each one of these channels are connected either to the periplasmic side of the membrane, 
having a higher electrochemical potential, or to the cytoplasm, characterized by a lower 
potential. The source (S) half-channel delivers protons from the periplasm (the positive, 
or, P-side of the membrane) to the rotor-stator interface, and, flnally, to the proton-binding 
sites in the rotor plane. The drain {D) channel translocates protons from this interface to 
the cytoplasmic, or negative, N-side of the membrane. The rotor sites can be protonated 
or deprotonated when they are in the immediate vicinity of the source or drain channels. 
In ATP synthase from Escherichia coli the residues cAsp-61 are equally placed along the 



rotating ring and serve as proton-binding sites carrying protons from the source to the drain 
channel llj. An empty rotor site has a negative charge. 



B. Motor mechanisms: Brownian ratchet and power stroke 

It is assumed [3, Q] that protons, dehvered to the left side (Fig. 1) of the rotor-stator 
interface by the source channel, bind and neutralize the empty rotor sites, and after that, 
travel clockwise (looking from the upper, periplasmic, side of the membrane, as shown in 
Fig. la). Notice that all rotor sites outside of the rotor-stator interface are occupied with 
protons and have no electric charge. The unprotonated negatively charged sites cannot face 
the lipid core of the membrane because of a large ( > 500 meV) desolvation penalty 17|] . The 
proton escapes from the rotor during the process of clockwise-directed rotor diffusion, when 
the corresponding site approaches the drain (cytoplasmic) channel located at the right side 
of the rotor-stator strip (Fig. 1). The negatively-charged empty site is not able to move back 
now because of the strong repulsion of the lipid medium. However, it can diffuse forward, in 
the clockwise (CW) direction, since there is no energy penalty for the charge placed inside 
the hydrophilic rotor-stator interface. The site is populated and neutralized again when it 
reaches the source channel at the left side of the interface, and the wheel can continue its 
CW rotation. 

Because of random thermal forces, the rotor can also rotate in the counterclockwise 
(CCW) direction, and the populated site would then approach the drain channel. In the 
most probable event, the site is depopulated and reflected from the lipid border of the 
rotor-stator interface. However, with a smaller probability, the neutral site, occupied with 
a proton, crosses the border, and the wheel rotates counterclockwise. When the population 
of the rotor sites occurs on the left side of the rotor-stator interface, and the depopulation 
happens on the right side, the rotor diffusion is biased to rotate CW. This means that the 

1^ 



Fq motor works as a Brownian ratchet 
References 



191, 



20|. 



2l| have shown, however, that the Brownian ratchet model is not able 
to explain the operation of the Fq motor in the presence of an external torque (about 41 
pN nm), which is produced by the Fi motor and acts in the opposite (CCW) direction. In 
addition to two half-channels, the stator contains a positive charge (a residue aArg-210 in 
E. coli) which can be located between two proton-conducting half-channels in the rotor plane 



[s], lljj. As schematically shown in Fig. la, upon passing the source channel, the rotor site 
releases its proton and becomes negatively charged. This charge is attracted by the positive 
residue, thereby generating a torque moving the rotor in the CW direction (power stroke 
component). Reference 8| has demonstrated that the positive stator charge is necessary 
for the efficient reverse operation of ATP synthase, when it works as a proton pump under 
the action of the external torque produced by ATP hydrolysis. The irreplaceable role of 
the aArg-210 residue for the ATP-driven pumping activity of the FqFi complex in E. coli 
has been proven by direct mutagenesis studies [22j. It should be noted that, despite the 
recent significant progress in understanding the operation of ATP synthase, the physical 
mechanism of the torque generation and proton translocation by the Fq motor calls for 
further investigation 12|, ll3| . Of special interest is the possibility to mimic the main features 
of this biological engine with the goal to create an artificial nanomachine with an almost 
perfect energy conversion. 



C. Quantitative modelling 



Here we exp 
rotary motors 



ore a simple model incorporating the most important aspects of biomolecular 
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111 and employing a quantitative treatment based on methods 



of condensed matter physics 



23 



24| . These approaches have been previously applied to 



nano-electromechanical systems (NEMS) with their mechanical motion affecting the electri- 
cal properties of electronic devices |25|. Similar processes take place in nanoscale biological 
objects, where electrical and mechanical degrees of freedom are also strongly coupled, mak- 
ing them living counterparts to artificial NEMS. Note that only nano-oscillators have been 
extensively studied by theorists, although a single-molecule rotor and nanoelectromechan- 
ical rotational actuators 26j have been demonstrated experimentally. To the best of our 



knowledge, no theoretical investigations of rotary NEMS have been reported yet. It can be 
expected that artificial nanoengines built on the operating principles of biological motors 
could achieve the same level of efficiency and performance. 

ATP synthase, along with many other enzymes {4], should be considered as a nonequilib- 
rium open system whose operation is supported by a permanent flow of protons. To describe 
this flow we attach the Fq motor to two proton reservoirs related to the periplasmic (positive, 
P) side of the membrane and to the cytoplasmic (negative, N) side of the membrane. In the 



case of the respiratory chain [l| , the P-side and the N-side correspond to the intermembrane 
space and to the matrix, respectively, divided by an inner mitochondrial membrane. Usually 
the proton transfer between the rotor sites and the P and N aqueous sides of the membrane 

transi- 



is represented by a set of Markov equations with a phenomenologically constructec 
jsl, isl, 21]. Here, with the methods of quantum transport theory jsl. H. 



tion matrix 
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25|, 



we derive the rate equations and provide simple expressions for transition coefficients, which 
explicitly depend on the difference between electrochemical potentials of the proton reser- 
voirs and on the distance between the proton-binding sites on the rotor and on the stator. 
This approach gives a clear physical picture of the loading and unloading events as well as 
the torque-generation process as a whole. 

The paper is organized as follows. In Sec. II we describe in detail the model and formulate 
the Hamiltonian of the system. In Sec. Ill we present a Langevin equation for the mechanical 
motion of the rotor and derive the rate equations for the populations of the proton-binding 
sites on the rotor. In Sec. IV we numerically solve the equations of motion and discuss our 
results. Section V is devoted to conclusions. 



II. MODEL 

We consider twelve equally-spaced proton-binding sites (cr = 1,...,12), attached to 
a ring-shaped rotor (c-subunit) of radius vq, at the points with the angular coordinates 
00- = cr ■ 00, where (po = n/Q is the angular distance between the rotor sites. Then, if the 
rotor is turned through an angle 0, the proton-binding site a has a coordinate 

r<7 = {ro sin(0 + 0^), cos(0 + 0^)}. 

We choose the coordinate origin at the center of the rotor wheel. As shown in Fig. 1, the 
stator part of the motor (or a-subunit) contains the source half-channel, S, which delivers 
protons from the P-side (periplasm) of the membrane to the point with an angular position 
05 = TT on the rotor-stator interface. The drain half-channel, D, also belonging to the stator, 
connects the point (j)£, = it — (pQ with the cytoplasmic (N) side of the membrane. Also, the 
stator has a positive charge q\e\ located in the rotor plane on the y-axis, with coordinates 
Tq = {0, — ro — d}. Here d is the distance from the charge q\e\ to the closest proton-binding 
site on the rotor, e is the electron charge. As a result, the Hamiltonian of the system has 



the form: 

ka 

E,n, + [f/g(0 + 0,) + f/con(0 + 0<x)] (1 - n,), (1) 

where the Hamihonian if^ governs the mechanical motion of the cyhnder-shaped rotator, 
characterized by the angle 0, counted clockwise from the y-axis. The occupation of the 
proton-binding site cr, having eigenenergy E^j, is described by the creation and annihilation 
operators of protons, a+,ao-, with a corresponding population operator n^j = a'^a^. Fermi 
operators c^^, Cka are related to the state, with momentum k, of the proton in the source and 
drain reservoirs {a = S,D) with energy Ej^a- The Hamiltonian -fftun describes the proton 
transfer between the proton-binding sites on the rotor and the source and drain channels on 
the stator. 



A. Coulomb interactions 

The empty site (with = 0) has a negative charge, and, because of this, it is attracted by 
the stator charge q\e\. The screened Coulomb coupling between the negative charge, — |e|, of 



the empty rotor site and the positive stator charge q\e\, separated 



distance rg(0 + (p^) = Ir^- — r^l, is determined by the potential 



jy the angular-dependent 



Uq{(^ + (^a) = - — — — exp , (2) 

47r£oerg(0 + 0^) V ) 

where Eq is the electric permittivity, e is the dielectric constant of the medium, is the 
Debye screening length, and 



rg(0 + 0a) = \jrl + (ro + df + 2ro(ro + d) cos((/. + 0^). (3) 

To take into account an energetic penalty for the charged rotor site cr, when it faces a 
lipid bilayer outside of the rotor-stator strip, we introduce a confinement potential, 

t^con(0 + 0<x) = f/c{l - exp[-Ae(l " COs(0 + 0, - 0c))']}, (4) 



where an angle 0c corresponds to the middle of the rotor-stator strip, 0^ = vr — 0o/2, a 
parameter Ac is inversely proportional to the width of the interface, and a potential U^. 



represents a desolvation penalty. The transfer of an ion with charge — |e| from water (with 
a dielectric constant £i) to the hydrophobic membrane (having a constant 62) is subject to 
the energy penalty |17l |. 



U, (meV) 



1440 e2 / 1 
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£1 



(5) 



2a \e2 

where a is the size of the cavity (in nm) over which the charge is spread. For ei = 80, 62 = 3, 
and a = 0.2 nm, the desolvation energy is about Uc = 1160 meV. We assume here that the 
proton-proton Coulomb interaction between sites is small enough, so that the loading and 
unloading of different sites occurs independently. 



B. Proton transfer 



Protons in the source and drain reservoirs are characterized by the Fermi distributions, 

-1 



/a(w) 



exp 



T 



+ 1 



(6) 



with temperature T {h= 1, = 1) and electrochemical potentials ns = V/2,fiD = — ^/2, 
where V is the proton voltage build-up, in units of energy, meV. For ATP synthase in Es- 
cherichia coli, protons can be translocated from the aqueous sides of the membrane to the 

n 

rotor-stator interface by a set of hydrogen-bonded chains [27| . The mechanism of a proton 
transfer between the rotor sites cAsp-61 and the terminal residues of the source channel, 
aAsn-214, and the drain channel, aSer-206, is not completely understood. However, molecu- 
lar dynamics simulations ll| demonstrate that hydrogen bonds can be formed between the 
proton-binding sites on the rotor (cAsp-61) and the terminal residues (aAsn-214 or aSer-206). 
These bonds are able to transfer protons within picoseconds either by collective tunneling 



271], or, which is more probable, by classical hopping 8|, [ill]. The proton translocation pro- 



cess strongly depends on the distance between the rotor sites and the stator residues and, at 
the molecular level, it can be facilitated by internal rotations of transmembrane helices as 



well as by the motion of other key elements UJ, l28|, |29l]. In our simplified model, we mimic 



the stator-rotor proton transitions with the Hamiltonian 



(7) 



kaa 



which describes the effective proton tunneling with amplitudes T^a, multiplied by factors 
w^a(0 + 0CT)- The factors Wa{(t> + (t>a) depend on the distance between the site a and the final 



residue of the a-channel located at the point = {rosiiK/)^, vq cos i^a}, 

\ra - = V2ro^/l - cos(0 + 0^ - 4>a). (8) 

We approximate this dependence by the exponential function, 

+ =exp(-Ar|r^-rc,|), (9) 

characterized by a steepness A,.. The value of is inversely proportional to the size of the 
molecular groups participating in the proton transport (either cAsp-61 and aAsn-214, or 
aSer-206). It should be noted that the specific functional form of the factor Waa{<P) is not 
crucial for the model under study. 



III. EQUATIONS OF MOTION 

The viscous medium creates a torque, 

X^C/^. (10) 

acting on the cylinder- shaped rotor, together with a stochastic force ^. For a cylinder with 
radius tq and height h, rotating in a medium with viscosity rj, the drag coefficient Cr is 



defined by the formula [30 



3l|: 



C. = 4n7]rlh. (11) 

It follows from the Hamiltonian ([T]), that the processes of loading and unloading of protons 
at the rotor-stator interface generates an additional torque. Thus, the following Langevin 
equation describes the biased overdamped diffusion of the rotor: 

= e + ^cxt - 5^ (1 - n^) [U,{(f) + + U,oni(f) + (f).)]- (12) 

(7 

Here ^(t) is a zero-mean value Gaussian fluctuation source, characterized by the correlation 
function: {^(t)^(t')) = 2T(rS(t—t'), and T^xt is an external torque produced by the Fi-motor. 
We expect that the Fq motor will rotate in the CW direction (looking from the periplasm side 
of the membrane). Therefore, the torque from Fi, which decelerates this motion, should have 
a negative sign, T^xt < 0. The effects of the rotor-stator proton transitions on the mechanical 
motion, resulting from the Hamiltonian Ht^n, are assumed to be negligibly small. 



The "chemical" part of the problem, namely, the process of loading and unloading the 
proton-binding sites on the rotor, is governed hy the Heisenberg equations for the population 
operators, = a^a^, which can be derived |5|, l2J] from Eqs. ([1]) and ([7]): 



ka 



where the reservoir operators, Cka, are represented as a sum of the free term, c[!^\ and the 
term describing the reservoir response: 

Cka = cfl - Tka'^ 1 dtigl.^{t,ti)a^{ti)waa{ti). (14) 

Here 

fkait^ti) = -iew[-iEka{t - ti)] 9{t - h) 

is the retarded Green function of the a-reservoir, 9(t — ti) is the Heaviside step function, 
and 

The correlator of the free reservoir operators is determined by the Fermi-distribution function 

fa{Eka)'- 

(ctl^HyfM) = UEka) exp[tEkait - ti)]. 

On substituting Eq. (fT^ into Eq. (fT3|) . and averaging the latter equation over the fluctuations 
of the proton reservoirs, we obtain 

5^|Tfc„|2 f dtiUEka)e'^''-^'-''^Waait)WaAh) X 

{[a,{t),at{h)\+)9{t-U) + h.c. (15) 
We introduce here the transition rates between the proton-binding sites and the reservoirs, 

k 

which are independent of energy in the wide-band limit approximation. Assuming weak 
coupling between the sites and the reservoirs, we employ free-evolving proton operators, 
a^{t) = e~*'^'^'^*~*^-'ao-(ti), to calculate the anticommutator in Eq. (ITSl) . Here 

= E,-t/,(0 + 0,)-[/,on(0 + 0.), (16) 



is the total eigenenergy of the proton on the site a including contributions of the stator 
charge potential, Uq{(f)), and the confining potential, Ucon{4>)- As a result, we derive a set of 
rate equations for the populations n^, partially averaged over the Fermi distributions of the 
proton reservoirs: 

ri<x + 5Zr,,(0) ri, = 5Zr,,(0) (17) 

a a 

Hereafter we drop the averaging brackets (...) and introduce the angular- dependent transi- 
tion coefficients, 

^aa{(f>) = wl{(f>(t) + (pa). 



The rate equations f[T7l) replace the phenomenological Markovian equations 9|, |l5|, |2l|] , which 
are usually employed for a description of the loading and unloading of the rotor sites. The 
characteristic time of proton transfer to and out of the proton-binding sites [s], [ll| is much 
shorter than the time scale of the rotation angle (j). Accordingly, we can average the stochastic 
Eq. (fT2|) over fluctuations of the proton reservoirs without averaging over the fluctuations of 
the mechanical heat bath. The partially averaged proton population n^, involved in Eq. ( |T2l) . 
depends on the local fluctuating value of the rotor angle 0(t). In the next section we solve 
numerically the stochastic equation (fT2!l together with the system of rate equations (fT7|l and 
investigate various regimes of the rotary nanomotor. 



IV. RESULTS 



We consider a cylinder-shaped motor with radius ro = 3 nm and height h = 6 nm, 
rotating in a medium with viscosity coefficient rj = 1 Pa-s, which is 1000 times higher 
than the viscosity of water. For a protein environment, surrounding the stator charge q\e\ 
and having a dielectric constant 6 = 3, the energy of the Coulomb coupling Uq{(j)) ([2]) is 
proportional to the factor qe"^ /{ATreos) = q x 480 meV-nm. The Debye screening radius is 
Ts = 1 nm. We also assume that both, the stator charge and the end of the source channel, 
are located in the rotor plane on the y-axis with = tt, and that the stator charge is 
offset by the minimum distance d = 0.9 nm from the circumference of the rotor ring. The 
drain channel is shifted to the right by the angle 0o : 0d = — 0o- For the steepness A,, 
involved in the function Wa{4>), we choose the value: A,. = (0.25 nm)~^. The confinement 
potential (jl]) is characterized by the parameters Uc = 1160 meV, Ac = 20. We assume 



that all proton-binding sites (cr = 1, . . . , 12) have the same eigenenergy, E„ = Eq. The 
proton transport through the system occurs if the 0-dependent energy levels of the proton- 
binding sites, E{(f)) = Eq — Ug{(f)) — Ucon{(p), located near the source and drain half- channels, 
= (pSi'pDi fit the transport window: 

^ > E{cPs) > Ei<PD) > - ^. (18) 

The energy Eq = — 90 meV gives the values E{(f)o) = — 84 meV, E{(f)s) = 100 meV, which 
meet these conditions at the transmembrane potential V = 250 meV. It is assumed that the 
time scales for transitions of protons between the reservoirs and the rotor sites are about 
= r^^= 0.5 /is. 

A. Rotations and site depopulations 

In Fig. 2a we plot a number of full rotations, 0(t)/27r, as a function of time (in ms), 
obtained from the numerical solution of the Langevin Eq. f|T2|) . coupled to the rate Eqs. f|T7j) . 
at the source-drain voltage V = 250 meV, the stator charge |e| (^'=1), and at temperature 
T = 300 K. We take into account here the constant load torque from the Fi motor, T^xt = ~ 
41 pN-nm, which is enough to produce three ATP molecules per rotation of the Fi motor. 

In the initial state all sites are occupied with protons, and, because of this, are electrically 
neutral. The process of site depopulation is illustrated in Fig. 2b, where we plot the 
probability to have the empty a-site, function of time for cr = 1, . . . , 12. 

The maximum value of n„ on this graph corresponds to the unoccupied proton-binding site 
a. The site "5" is depopulated in the first turn, because it is initially in the closest position 
to the drain channel. The empty site is attracted by the positive stator charge, and this 
force turns the rotor an angle tc/Q (see Fig. 2a). The site "5" is occupied again when it 
reaches the source channel. In Fig. 2c we plot the number Ns of protons, transferred from 
the source reservoir, as a function of time. It is evident from Fig. 2c that loading the site 
"5" is accompanied by the transfer of a proton from the source channel. The depopulation 
of the site "4" begins at the same time as the depopulation of the site "5". This process 
repeats over and over, resulting in a continuous unidirectional rotation of the rotary ring. 
The loaded motor makes a complete turn in a time interval ~ 0.25 ms, which corresponds to 
the frequency of rotations near 4 kHz; in so doing the system translocates about 12 protons 



down the potential gradient of 250 meV. 

It should be noted that here we have an asymmetric configuration, where the energy level 
of the proton near the source channel is higher than the proton energy near the drain outlet, 
E{(ps) > E{(f)£)). The energy drop, 

AE = E{<Ps) - E{<Pd) = Ugi<PD)-U,i<Ps), (19) 

is directly converted into mechanical energy by the Fq motor. We note that Uq{(j)) < 
and f/con(0s) = ^con(0D)- At the opposite sign of the transmembrane potential difference, 
V < 0, the motor does not produce any unidirectional rotation, since in this case the proton 
energies do not fit the transport window (ITSl) . 

B. Torque generation and shuttling 

In Fig. 3 the average rotational frequency, {Q)/2'k = {(j)) /27r, shown in (a), and the proton 
current. Is = Ns/tr-, (b), averaged over the time interval r/j = 2 ms as well as over several 
realizations, are presented as functions of the eigenenergy of the proton-binding states on 
the rotor, Eq, at three values of the temperature, T = 4.2 K (liquid ^He); T = 77.2 K (liquid 
^^N); and T = 300 K. The other parameters, namely, the proton voltage, V = 250 mV, the 
normalized stator charge, q = 1, and the external torque, T^xt = — 41 pN-nm, correspond 
to the loaded Fq motor. At resonant values of Eq : Eq — 90 meV, and at high enough 
temperatures, T > 50 K, the motor works in a torque-generating regime characterized by 
a positive (CW) direction of rotation, (Q) > 0, and positive current. Is > 0, when protons 
fiow downhill, from the source to the drain reservoir. In this regime both mechanisms, the 
power stroke and the Brownian ratchet, contribute to the torque generation. 

At low temperatures, e.g. when T = 4.2 K, the ratchet mechanism practically turns 
off, and the system is not able to produce enough torque to execute the full-circle rotation 
in the presence of the counteracting load torque T^xt- In this shuttling mode, one proton- 
binding site oscillates back and forth between the source and drain channels, translocating 
the protons {Is > 0), but no unidirectional rotation is generated (Q) = 0). 

When Eq < — 160 meV the system is out of the transport windows, and the proton- 
binding sites on the rotor are always populated. Then, the rotor follows the negative external 
torque and rotates in the CCW direction without transferring protons (i.e., (Q) < 0, Is — 0). 



C. Torque generation and proton pumping 



In Fig. 4 we show the dependence of (a) the average speed of rotations, {Q)/2'7i, and (b) the 
average particle current, Is (average number of protons transferred from the source channel 
per one millisecond), on the voltage V for the whole range of the CCW-directed external 
torque T^^t, from zero to — 120 pN-nm, at Eq = — 90 meV, and T = 300 K. For values 
of the external torque between zero and — 80 pN-nm and at a sufficiently high voltage, 
V > 220 meV, the motor performs a CW rotation (Fig. 4a) with a maximum frequency 
~10 kHz (the torque- generating mode). In the process, protons flow downhill, from the 
source (/i5 = V/2) to the drain channel (/iz? = ~V/2). The CCW- directed external torque 
hampers this motion, and when |%xtj> 80 pN-nm the rotation stops (Fig. 4a). 



In the regime of ATP hydrolysis [13|], the Fi motor produces a sufficient torque to drive 
the rotation of the Fq motor in the reverse (CCW) direction. In our case this regime takes 
place when |%xt| > 50 pN-nm, provided that the transmembrane potential is small enough, 
V < 170 meV (Fig. 4a). It is evident from Fig. 4b that in this range of parameters the 
system works as a proton pump {Is < 0), which translocates protons against the gradient 
of the electrochemical potential with a pumping efficiency 

Effpump = ^^f^ - 20% (20) 

when V = 100 meV and T^xt = — 108 pN-nm. 

Note, that the voltage threshold for the torque-generating regime {V > 210 meV when 
T^xt = — 41 pN-nm) corresponding to the value of the electrochemical potential difference 
across the bacterial membrane and the inner mitochondrial membrane, is not preassigned 
a priori, but it is naturally determined by the configuration parameters of the model; for 
example, by the charge and the position of the stator as well as by the value of the load 
torque necessary to drive the ATP synthesis. 

D. Contributions from the Brownian ratchet and power stroke components 

The model studied here can be considered as a combination of Brownian ratchet and 
power stroke components 10, l2l|]. The Brownian ratchet harnesses the energy of thermal 
fluctuations, biased by chemical reactions, to generate a unidirectional rotation. In our case, 
the ratchet component is introduced via a fluctuation force ^, working in combination with 



the processes of loading and unloading protons in the presence of the confining potential 
f^on (see Eqs. f|T2l) . f|T7|) ). It is expected that the ratchet contribution should be significantly 
diminished at low temperatures. 

The power stroke component is due to the Coulomb attraction between the positive stator 
charge, g|e|, and the negatively charged unoccupied rotor site. This part of the torque- 
generating process is proportional to the coefficient g, which is the stator charge, measured 
in units of |e|. For the torque-generating and pumping process in the Fq motor of E. coli, 
the importance of the residue cArg-210, carrying a positive charge, has been emphasized in 
Refs. y. 
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In Fig. 5, taking the values d = 0.9 nm, Eq = — 90 meV, and V = 250 meV, we plot: (a) 
;he average torque (T) generated by the isolated motor (at |7^xt| = 0), and (b) the efficiency 

3, 

Eff ^ (21) 

of the system as a function of the stator charge q, normalized by |e|, at low (T = 4.2 K, 
dashed blue line), intermediate (T = 77.2 K, dash-dotted green line), and high (T = 300 K, 
continuous red line) temperatures. 

In the pure Brownian ratchet regime (no stator charge, q = 0) the system generates a 
pronounced torque, (T) ~ 12 pN-nm, if the temperature is high, T = 300 K. However, this 
torque is not enough to overcome the load torque, T^xt = — 41 pN-nm, which is necessary 
for ATP synthesis. For instance, at liquid Helium temperatures (T = 4.2 K) the Brownian 
ratchet component is strongly suppressed. 

Nevertheless, the power stroke mechanism can generate the torque (T)inax — 60 pN-nm, 
which is higher than the load torque from the Fi motor. The isolated Fq motor demonstrates 
this peak value of the torque and the efficiency, Effmax — 80%, at the stator charge near |e| 
(g^l). 

It follows from Fig. 3 that the power stroke component alone is not able to drive ATP 
synthesis, because at the load conditions, T^xt = — 41 pN-nm, and at T = 4.2 K, the 
system operates in the shuttling mode, with a proton current but generating no unidirectional 
motion. Notice also that despite the presence of a medium with high viscosity coefficient, 
rj = 1 Pa-s = 1000 ■ rj^ater, and despite the very small size, tq = 3 nm, the motor performs 
quite well and generates a significant torque at physiologically reasonable parameters. 

Note that, because of the asymmetric configuration, this motor does not produce a unidi- 



rectional torque for a negative electrochemical potential difference {V < 0), but it switches 
the direction of rotation if the drain channel (see Fig. 1) is placed on the other side of the 
source channel (at (pn = 7r + 0o)- This feature, as well as the nanoelectromechanical method 
as a whole, can be useful for explaining the switching ability of bacterial flagellar motors, 
which allows the chemotaxis of bacteria jl^ . . 



V. CONCLUSIONS 



We have examined a simple nanoelectromechanical model of a proton-driven rotary 
nanomotor which mimics the basic operating principles and the configuration of the Fq- 
motor of ATP synthase in bacterial and mitochondrial membranes. Treating ATP synthase 
as a nonequilibrium open system, coupled to two proton reservoirs, we have derived a set 
of rate equations, which describes the loading and unloading of the proton-binding sites on 
the rotor portion of the motor. 

At normal conditions, the isolated motor generates a torque of about 60 pN-nm with 
efficiency near 80%, whereas the motor working against a constant load, created by the Fi 
portion, exhibits a unidirectional rotation (Fig. 2) with frequency ~ 4 kHz. It is shown 
that, depending on the temperature T, the external torque T^xt, the proton voltage build-up 
V, and the energy Eq of the rotor sites, the system operates in three different regimes: the 
torque-generating mode; the shuttling regime; and the pumping mode (see Fig. 3 and Fig. 4). 

In the torque-generating mode, energy (stored in the gradient of the proton electrochem- 
ical potential) is converted into mechanical energy of the rotor. This energy transfer drives 
the unidirectional rotation of the motor, (Q) > 0, and the positive particle current. Is > 0, 
corresponding to the downhill flow of protons from the periplasmic (P) to the cytoplasmic 
(N) side of the membrane. 

In the shuttling regime, which takes place at low temperatures, the rotor vibrates with 
a small amplitude near the initial point without performing a full circle. In doing so, 
the proton-binding site, located between the source and drain half-channels, works as a 
nanomechanical shuttle carrying protons from the P-side to the N-side of the membrane. 
This regime is distinguished by zero-frequency rotations, (Q) = 0, and by a positive proton 
current, Jc > 0. 

In the „ .node Q. ATP synthase operates in reverse, where ATP hydrolysis 



drives rotations of the Fi motor. This external torque, %xt; is transmitted to the Fq motor, 
which pumps protons uphill, from the N to P membrane side. In our case, the CCW-directed 
external torque, T^xt = — 120 pN-nm, is enough to translocate, in one millisecond, about 90 
protons against the electrochemical gradient of 100 meV (see Fig. 4b). 

We have studied the performance of the model at different temperatures and found that 
this performance depends significantly on the Coulomb interaction, Ug, between the positive 
stator charge and the negatively charged empty rotor site. The torque produced by the 
Brownian ratchet mechanism alone is not sufficient to overcome the load torque from the 
Fi motor. A key feature of the model is that, because of Coulomb coupling to the stator 
charge, the energy of the proton on the rotor site, which is close to the source half-channel, 
E{(f)s), is higher, than the proton energy of the site located near the outlet of the drain 
channel, E {(()£,). 

For positive proton voltage V , when the electrochemical potential of the source, = ^/2, 
exceeds the potential of the drain, fio = —V/2, the rotor sites can be loaded with protons 
at the source and unloaded at the drain, if the energies E{(f)s) and E{(j)u) fit the transport 
window: fis > E{(f)s) > -E(0d) > /^d- This energy difference, AE = E{(f)s) — E{(f)s) = 
Uq{4>D) — Uq{(ps), is directly converted into mechanical energy. It is shown, however, that 
at low temperatures, when the Brownian ratchet component is suppressed, the motor works 
in the shuttling mode, and no torque is generated. This means that, in agreement with 
previous results fl, Q, the power stroke and the Brownian ratchet components should 
work together for the efficient operation of the system. These conclusions can be applied 
both to biological rotary motors and to artificial nanoengines. Notice also that an efficient 
and powerful future synthetic motor based on this bio-inspired design would be expected 
to withstand a more severe environment, in particular, much lower temperatures, than 
its biological counterparts. The present design can be compared with modern industrial 
electrical motors having an efficiency about 90%. However, it should be emphasized that 
this efficiency decreases drastically when decreasing the size of the motors. Nano-scale 
motors operating in a warm and very viscous medium, even with an efficiency 80%, are 
worth building and studying as a function of various operating conditions. 
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FIG. 1: (Color online) Schematic diagram of a rotary nanomotor powered by an electrochemical 
potential gradient, (a) shows the top view from the periplasmic (P) side of the membrane; (b) 
shows the side view. Panel (a) shows the rotor, which has twelve sites which rotate as a Ferris 
wheel. In the example shown in (a), the site number five provides a proton to the drain (D) 
channel of the stator, and the proton then moves to the cytoplasm, as shown in (b). This now 
negatively-charged site number five in (a) feels a Coulomb attraction from the positive (red) stator 
charge, which propels the wheel to rotate an angle 0o = 30°. After site number five moves to the 
location of the next site, the source (S) donates a proton to that site, neutralizing it, and the 30° 
rotation process is ready to start all over again. In (b), the protons, indicated by "+", move from 
the periplasm (P), which acts as a source (S), to the cytoplasm {N), acting as a proton drain (D), 
on the other side of the membrane. In (a) and (b), an orange arrow represents the Coulomb force 
between the empty rotor site and the positive stator charge. 




FIG. 2: (Color online) (a) Time dependence of the number of full rotations (j){t)/2TT at ^ = 250 
meV, T^xt = — 41 pN-nm, Eq = — meV, and at T = 300 K; (b) Depopulations of the proton- 
binding sites, no- = (1 — no-), o" = 1, . . . , 12, versus time (in ms); (c) The number of protons, Ns, 
transferred from the source to the drain reservoir, versus time (in ms). Notice the periodicity in 
(b). It follows from Fig. 2a, that the motor performs a bit more than one full CW rotation for 
a period of time ~ 0.25 ms. This rotation is accompanied by the cyclic loading and unloading of 
the proton-binding sites on the rotor (Fig. 2b), starting with the site number five (see Figs, la and 
2b). All sites were populated at the initial time t = 0. The rotor works as a Ferris wheel carrying 
protons from the source reservoir S to the drain D. The number of protons, A^s, translocated form 
the source to the drain, increases stepwise during the rotation of the wheel (Fig. 2c). Almost 12 
protons are transferred through the 12-site system after the full circle, as seen in Fig. 2c. 




FIG. 3: (Color online) (a) The average frequency of rotations, (0) jlix (in kHz), and (b) the average 
proton current Is (the number of protons, transferred through the system in one millisecond) as 
functions of the energy of the rotor sites, Eq (in meV), at F = 250 meV, T^xt = — 41 pN-nm, 
and at different temperatures: T = 4.2 K; T = 77.2 K; and T = 300 K. In the torque- generating 
regime (T = 300 K, i?o ~ —80 meV), despite the counteracting load torque T^xt, the system 
rotates CW, (17) > 0, performing about 5 full circles (Fig. 3a, red continuous line) and carrying 
more than 60 protons from the source to the drain in one millisecond (Fig. 3b). However, at very 
low temperatures, T = 4.2 K, the Brownian ratchet component of the torque-generating process 
turns off, and the rotor is not able to overcome the external torque. The system is now in the 
shuttling regime oscillating near a fixed position with zero average frequency of rotations, (O) = 0, 
but carrying about 20 protons per millisecond (Figs. 3a and 3b, flat plateaus in the blue dashed 
curves) . 
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FIG. 4: (Color online) (a) the average speed of rotations (J7)/27r (in kHz), and (b) the average 
proton current, Is, as functions of the proton voltage, V (in meV), and the external torque, T^xt (in 
pN-nm), at £"0 = — 90 meV, and T = 300 K. Figure 4a shows that for a low enough counteracting 
load torque, |7^xt| < 80 pN-nm, and at a high proton voltage, V > 220 meV, the motor rotates 
clockwise with a maximum frequency of about 10 kHz. This rotation is accompanied by the positive 
proton current reaching the value Is ~ 100 protons/ms (Fig. 4b). The blue (negative) regions in 
Fig. 4 correspond to the regime of ATP hydrolysis (pumping mode), when at low voltages, V < 170 
meV, the external torque, |7^xt| > 50 pN-nm, is powerful enough to drive the motor rotation in 
the CCW direction ((i7) < 0), pumping protons (Is < 0) from the drain to the source against 
the proton electrochemical gradient V (Fig. 4b). 





FIG. 5: (Color online) (a) Average torque, (T) (in pN-nm), and efficiency, Eff, of the isolated motor 
{%ixt = 0) as functions of the normalized stator charge q (in units of |e|) for V = 250 meV, Eq = 
— 90 meV, and for three different temperatures: T = 4.2 K (blue dashed line); T = 77.2 K (green 
dash-dotted line); and T = 300 K (red continuous line). Without the stator charge (Fig. 5a, red 
continuous line, T = 300 K, g ~ 0) the ratchet mechanism produces the torque ~ 12 pN-nm, which 
is not enough to drive ATP synthesis. Near the optimal value of the stator charge {q ~ 1), but 
at very low temperatures, T = 4.2 K, the isolated motor generates the torque 60 pN-nm by means 
of the power stroke mechanism (Fig. 5a, blue dashed line). This torque is higher than the CCW 
torque (— 41 pN-nm) necessary for ATP production. However, at these conditions, the Fq motor 
coupled to the load (T^xt = — 41 pN-nm) switches from the torque-generating mode to the shuttling 
regime, as shown in Fig. 3, and no average torque is produced. Figure 5b demonstrates that the 
efficiency of the motor can reach ~ 80%, when both components contribute to the torque-generating 
process. 



